The role of Zhang-Rice singlet-hke excitations in one-dimensional cuprates 
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We present the first calculation of the electron-energy loss spectrum of infinite one-dimensional 
undoped CuOa chains within a multi-band Hubbard model. The results show good agreement with 
experimental spectra of Sr2Cu03. The main feature in the spectra is found to be due to the formation 
of Zhang-Rice singlet-like excitations. The q-dependence of these excitations is a consequence of 
the inner structure of the Zhang- Rice singlet. This makes the inclusion of the oxygen degrees of 
freedom essential for the description of the relevant excitations. We observe that no enhanced 
intersite Coulomb repulsion is necessary to explain the experimental data. 

PACS numbers: 71.27.-fa, 71.45. Gm, 71.10.Fd 
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Recently, charge excitations in the quasi one- 
dimensional compound Sr2Cu03 have been investigated 
both experimentally^"'^ and theoretically.^"^ Sr2Cu03 is 
composed of chains formed by CUO4 plaquettes which 
share the corner oxygens. The magnetic properties of 
these chains have been successfully described using a one- 
dimensional spin-i Heisenberg antiferromagnet.^"^ 

Experimentally, the electron-energy loss spectrum 
(EELS) of Sr2Cu03El shows several interesting features 
(see Fig. |l|): For small momentum transfer {q — 0.08 
A"-'^) parallel to the chain direction, one observes a broad 
peak around 2.4 eV energy loss, and two relatively sharp, 
smaller maxima at 4.5 and 5.2 eV. With increasing mo- 
mentum transfer, the lowest-energy feature shifts towards 
higher energy, reaching 3.2 eV at the zone boundary 
{q = 0.8 A~^). Thereby its spectral width decreases. 
In addition, the peaks at 4.5 and 5.2 eV lose spectral 
weight as the momentum transfer increases, while some 
less well-defined structures emerge around 6 eV. 

So far, these results have been compared only to cal- 
culations in an extended one-band Hubbard model. '^'^ 
From this comparison, Neudert et al.^ concluded that in 
Sr2Cu03 there is an unusually strong intersite Coulomb 
repulsion V: In the one-band model it is of the order of 1 
eV. It is argued that this large value of V allows for the 
formation of excitonic states which are observed in the 
experiment. One of the aims of this paper is to show that 
no intersite Coulomb repulsion is necessary to explain the 
basic features of the experiment, if the O degrees of free- 
dom are taken into account within the framework of a 
multi-band Hubbard model. 

We investigate the EELS spectrum of a one- 
dimensional CuOs chain system, using a multi-band 
Hubbard Hamiltonian at half-filling. In the hole picture 
this Hamiltonian reads 
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where c?|^ (Pja) create a hole with spin a in the i-th Cu 
3d orbital {j-th O 2p orbital), while nf^ ('^^o-) 
corresponding number operators. The first and second 
term in Eq. (nl) represent the atomic part of the Hamilto- 
nian, with the charge-transfer energy A, and the on-site 
Coulomb repulsion Ud between Cu 3d holes. The last two 
terms in Eq. (^ are the hybridization of Cu 'id and O 
2p orbitals (hopping strength tpd) and of O 2p orbitals 
(hopping strength tpp). The factors and give 
the correct sign for the hopping processes. Finally, (ij) 
denotes the summation over nearest neighbor pairs. 

The loss function in EELS experiments is directly pro- 
portional to the dynamical density-density correlation 
function Xp(ci;,q),llj which depends on the energy loss 
Lu and momentum transfer q. Xp('^i q) is calculated from 



Xp(c.,q) = - / dt e-^-'(vI/|[p_q(0),Pq(t)]|vI/) , (2) 
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where l^*) is the ground state of H, and pq is the Fourier 
transformed hole density. The ground state \^) is ap- 
proximated as follows :Lj We start from a Neel-ordered 
state I^at) with singly occupied Cu 3d orbitals (with al- 
ternating spin direction) and empty O 2p orbitals. Fluc- 
tuations are added to using an exponential form 



\^) = exp ^E ^o^F^.y^ 



(4) 



The fluctuation operators Fi^a describe various delocal- 
ization processes of a hole initially located in the Cu 3c? 
orbital at site i, wheipe-a summation over equivalent fi- 
nal sites takes placeO The parameters in Eq. (^) 
describe the strength of the delocalization processes and 
are determined self-consistently by solving the system of 
equations = 0, where L is the Liouville op- 

erator, defined as LA — \H, A\ for any operator A. These 
equations have to hold if \^) is the ground state of H. 
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Using Eqs. (||) and (Q), we calculate the EELS spefe. 
trum by means of Mori-Zwanzig projection technique.t£l 
For a set of operators D^, the so-called dynamical vari- 
ables, the following matrix equation approximately holds 
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'MDlj—^D,\^) = mDlD^l^) , (5) 



where z ~ uj+iO. In Eq. (|5|) the set of dynamical variables 
was assumed to be sufficiently large so that self-energy 
contributions can be neglected. The set {D^} contains 
the dynamical variable Dq = pq. Therefore, by solving 
Eq. an approximation for Eq. (H) can be obtained. 
Besides Dq, the set includes Da = Pq-Fo,Q for all a. The 
Fo^a are the fluctuation operators used in the ground 
state Eq. (Q), without the summation over equivalent fi- 
nal sites. We use altogether 12 dynamical variables and 
observe good convergence of the spectral function. 




FIG. 2. The theoretical results for the dominant exci- 
tation at 2.4 — 3.1 eV with a broadening of 0.02 eV. The 
momentum dependence of the spectrum is due to two dif- 
ferent effects. First, with increasing q the spectral weight 
shifts from excitation (a) to (b). Second, the energies are 
q-dependent. Both effects contribute to about one half of the 
full momentum dependence. 
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FIG. 1. Comparison of experimental data for Sr2Cu03 
(left), taken from Ref. 3, and the present theoretical results for 
the one-dimensional multi-band Hubbard model (right) . The 
theoretical line spectra have been convoluted with a Gaussian 
function of width 0.1 eV. For details see the text. 

In Fig. |l] the obtained results are compared to the ex- 
perimental spectra from Ref.0. The parameters in the 
Hamiltonian are chosen as follows: Ud = 8.8 aV, and 
tpp — 0.65 eV are kept constant at typical values.ll3 The 
values of A = 4.3 eV and tpd = 1.5 eV have been adjusted 
to obtain the correct position of the lowest energy feature 
at 2.5 eV for q = 0.01 A^^, and at 3.1 eV for q = 0.7 A'^. 
Thus, we effectively use only two free parameters. It is 
found that the value of A dominates the excitation en- 
ergy, which increases with increasing A. The dispersion 
of the peak depends mainly on tpd with increasing dis- 
persion for increasing hoppinEJ)arameter. As compared 
to the standard value 1.3 eVjlj tpd — 1.5 eV is slightly 
enhanced, in agroement with recent results of band struc- 
ture calculations.EJ 



The theoretical spectra consist of two excitations. The 
dominant excitation is at 2.45 eV for q = 0.1 A~^, and 
shifts to 3.05 eV for q — 0.7 A~^. Besides, a second 
excitation appears at 6.4 eV which has no dispersion. 

The low energy peak structure is shown in more detail 
in Fig. 2, where a smaller peak broadening has been used. 
As will be explained belowjjnainly two different Zhang- 
Rice singlet-like excitationsEJ lead to this peak structure. 
The q-dependence of the spectrum is due to two effects. 
Firstly, one observes a shift of spectral weight with in- 
creasing q between two excitations labelled with (a) and 
(b) in Fig. 2. Secondly, with increasing q the energies of 
the two peaks shift to higher values. 

The shift of spectral weight can be attributed to dif- 
ferent delocalization properties of the two final states. 
The excited state (a) in Fig. 2 which dominates the spec- 
trum for small momentum transfer is rather extended, 
see Fig. 3(a). This state has a rather small probability 
for the hole at its original plaquette. With increasing q 
the spectral weight shifts to another excited state, shown 
in Fig. 3(b), with a higher probability for the hole on its 
original Cu-site. This means that the character of the 
excitation changes from an extended to a more localized 
one, while still forming a Zhang-Rice singlet. 

This behavior can be understood by analyzing the rel- 
evant expectation values in Eq. (5) . For small values of q 
the frequency term {'i'\Djj^CDi,\'^) can be approximated 
by expanding e^'^'" « 1 + iqr in Eq. (||). This gives 

{'^\^o.h^^Om\'^) X q(r^ — Tj,) which is proportional to 
the fluctuation distance, thus favoring far-reaching exci- 
tations. This picture changes for large values of q, where 
stronger oscillations of the phase factor lead to a cance- 
lation of extended excitations. The result is a transfer of 
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spectral weight from delocalized towards more localized 
excitations with increasing q. 




FIG. 3. Hole delocalization properties of different excited 
states. Larger (smaller) dots symbolize a larger (smaller) den- 
sity of the hole originally located at the central plaquette. 
States (a) and (b) are Zhang- Rice singlet-like excitations with 
different delocalization properties. Excitation (a) has the 
largest spectral weight for small q, whereas excitation (b) has 
dominant spectral weight for large q, see also Fig. 2. Part (c) 
shows the local excitation at 6.4 eV, where the hole surrounds 
the central Cu site. 

The q-dependence of the energies, on the other hand, 
is a consequence of the inner structure of the Zhang-Rice 
singlet-like excitations. In both excitations (a) and (b) 
a hole hops onto the Cu site of its nearest neighbor pla- 
quette, see Fig. 3. Due to the Coulomb repulsion Ud, the 
hole which had originally occupied this Cu site is pushed 
away onto the surrounding O sites. Depending on the 
direction of this delocalization, this process leads to a 
q-dependence of the excitation energy. _ 

Next, we want to stress that the claim inl3 for the one- 
band Hubbard model that only the inclusion of the next- 
neighbor repulsion leads to the possibility of the forma- 
tion of an excitonic state is not consistent with our re- 
sults. In the one-band model such a repulsion leads to a 
binding energy of empty and doubly occupied sites due 
to the reduction of neighboring interactions. This bind- 
ing energy is proportional to V . However, as can be seen 
from ci^ct diagonalization calculations in the one-band 
model,ll3 the intersite repulsion mainly leads to an en- 
ergy shift of the EELS spectra. Thus, the parameter V 
in the one-band model serves only to adjust the energet- 
ical position of the spectra, and is not necessary in more 
realistic models. In the multi-band model, the formation 
of an exciton is only driven by the energetically favored 
formation of a Zhang-Rice singlet, and no further inclu- 
sion of next-neighbor repulsion is necessary. 



The important role of the Zhang-Rice singlet formation 
has been studied previously |a^o in an effective model for 
excitons in the Cu02 plane.llZl Like the one-band model, 
this effective model neglects inner degrees of freedom of 
the Zhang-Rice singlet. If this model is reduced to the 
CuOa chain, q-dependent energies are only possible for 
a non-vanishing O on-site Coulomb repulsion Up ^ 0. In 
contrast to these results, we find q-dependent energies 
for Up = 0. As described above this effect cannot be ex- 
plained in a model which neglects the inner structure of 
the singlet. 

Thus, our results show that both an inclusion of the 
0-sites and a complete description of the excitation is 
necessary to obtain the full dispersion. The 0-sites are 
essential for the correct description of the different char- 
acters of the singlet excitations, which leads to the shift 
of spectral weight from one excitation to another. On the 
other hand, taking account of the inner degrees of free- 
dom of the Zhang- Rice singlet leads to the q-dependence 
of the energies. 

The results of the projection technique do not correctly 
describe the experimentally observed width of the peak 
for small momentum transfer. A possible explanation 
is that not all excitations are included in the projec- 
tion space. The above discussion suggests that the width 
should be due to the presence of additional delocalized 
excitations. Processes which are neglected in the present 
calculation involve less important multiple excitations of 
holes beyond their original plaquette. 

Finally, although they are not the focus of this paper, 
we discuss some high-energy features. The excitation at 
6.4 eV in the theoretical spectra is due to a local pro- 
cess on the plaquette itself. Here, the hole is excited 
to the surrounding O sites, without leaving its original 
plaquette, see Fig. 3(c). The energy of this structure 
does not shift as a function of momentum transfer. Once 
again, a transfer of spectral weight towards this local- 
ized excitation with increasing values of q is observed. 
The plaquette excitation has a highly local character. 
Therefore, its spectral weight increases as a function of 
q compared to the more delocalized Zhang-Rice singlet 
excitations. For small q the spectral weight of the pla- 
quette peak is about 6 times smaller than that of the 
Zhang-Rice peak. As q increases, this ratio increases to 
about one half. One should note that the experimental 
spectra show no obvious features above 6 eV. However, 
since many different orbitals may contribute in this en- 
ergy range, we cannot expect a realistic description using 
a model that contains only Cu 3d and O 2p orbitals. This 
applies also to the experimental structure around 4.5 eV 
for small momentum transfer, which is not described by 
the present model. We assume that this feature is due to 
excitations which involve Sr orbitals, as has been argued 
before E 

In comparison with earlier works on Cu2p3/2 X-ray 
photoemission spectroscopy usiirg the same theoretical 
approachjEa we find that the character of the excitations 
in both experiments is very similar. Zhang-Rice singlet 
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and local excitations play an important role. In both 
experiments the dominant excitation at low energies is 
associated with a Zhang-Rice singlet formation. 

In conclusion, wc have carried ont the first calcula- 
tion of the EELS-spectrum for the one-dimensional CuOa 
chain by using a multi-band-Hubbard model. Our re- 
sults arc in good agreement with experimental results 
for Sr2Cu03. We find that the main feature in the spec- 
tra is due to the formation of Zhang-Rice singlet-like ex- 
citations. The momentum dependence of the spectrum 
is due to two effects. First, there is a shift of spectral 
weight from less localized to more localized final states. 
Second, the excitation energies are q-dependent. This 
q-dependence is found to be a consequence of the inner 
structure of the Zhang-Rice singlet. Therefore, the in- 
clusion of the O degrees of freedom is essential for the 
description of the relevant excitations. This has two im- 
portant consequences. Firstly, only a multi-band model 
allows the correct description of charge excitations. And, 
secondly, if a multi-band model is used, no intersite 
Coulomb repulsion is necessary. Furthermore, we observe 
the existence of a local excitation at large q- values. 
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